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Abstract. In this paper we show preliminary results concerning the 
chemical and energetic enrichment of the ICM by means of supernova- 
driven wind models in elliptical galaxies. These models are obtained by 
taking into account new prescriptions about supernova remnant evolution 
in the interstellar medium. We find that models, which can reproduce 
the Fe abundance and the [a/Fe] ratios observed in the ICM, predict 
that the SN energy input can provide about 0.3 keV per ICM particle. 
We have obtained this result by assuming that each SN explosion inject 
on the average into the ISM no more than 20% of its initial blast wave 
energy. The predicted energy per particle is not enough to break the 
cluster self-similarity but is more than predicted in previous models. 



1. Introduction 

Hydrodynamical simulations and semi-analytical models for clusters of galaxies 
need an extra-energy of about 1 keV per ICM particle (e.g. Wu et al. 2000; 
Bower et al. 2000; Bialek et al. 2000; Borgani et al. 2001) to be reconciled with 
X-ray observations. In this framework one of the open problems is to assess 
whether SNe can provide this energy in order to break the cluster self-similarity 
(e.g. Finoguenov et al. 2001) or not (e.g. Valageas & Silk 1999; Bower et al. 
2000). On the other hand, as shown by Matteucci (this conference), several 
successfull models for the chemical enrichment of the ICM have been developed 
in last twenty years. We remind that a good model for the ICM chemical enrich- 
ment must produce (Fe/H)/(Fe/H) Q ~ 0.3 (e.g. White 2000) and [a/Fe] ~ 
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(e.g. Ishimaru & Arimoto 1997; Renzini 2000). Since the galactic ejecta carry 
out energy with them we adopted the approach by Matteucci & Vettolani (1988) 
for computing the wind energy in models which match the observational con- 
straints on the chemistry of the ICM. We modified the chemical evolution codes 
developed by Matteucci & Gibson (1995) and Martinelli et al. (2000) includ- 
ing a new cooling time for supernova remnants (SNR) depending on metallicity 
(Cioffi & Shull 1991) as well as new prescriptions for SNIa cooling (Recchi et al. 
2001, section 2). Then we followed the evolution of the ICM abundances with 
redshift (section 3). In section 4 some results are shown and some conclusions 
are drawn. 



2. Chemical Evolution Model 

Chemical evolution models (both one and multizone) are those of Matteucci k, 
Gibson (1995) and Martinelli et al. (1998) . We focus here on the new energetic 
prescriptions implemented in the chemical evolution code. We follow the SNR 
evolution in the interstellar medium according to the cooling time by Cioffi & 
Shull (1991): 

t cooi = 1.49a0 4 efV 4/7 C 5/1 V, (1) 

where Q = Z/Z & , uq is the hydrogen number density, eo is the energy released 
during a SN explosion in units of 10 51 erg and we take always eo = 1- Old (Cox 
1972, long dashed & dotted line) and new (Cioffi & Shull 1991) cooling times are 
compared in figure 1, where metallicity and density of the ISM evolve in a self- 
consistent way as a functions of time. The new cooling time (solid line) is about 
3 times lower than the older one after 0.1 Gyr from the beginning of galactic 
evolution and that, soon after 0.2 Gyr, the metallicity becomes oversolar and 
makes the cooling ~ 10 times more efficient. The most important consequence 
is that, while in previous works flat initial mass functions (IMF) were preferred 
to get the right amount of elements in the ICM, now we show that models with 
Salpeter IMF as the best ones. The reason is that the faster the metallicity 
grows the more efficient the cooling is. So, if we compare a galaxy with a flat 
IMF (x=0.95) with a galaxy of the same initial mass with Salpeter IMF, the 
latter undergoes galactic wind earlier than the former and consequently eject 
more metals into the ICM. 

The second fundamental hypothesis on energetics of the ISM is that SNe la 
are allowed to transfer all of their initial blast wave energy into the ISM: 

tSNIa = 1 • (2) 

In fact Recchi et al. (2001) showed that radiative losses for SNIa are negligi- 
ble, since their explosions occur in a medium already heated by SNII. Therefore 
for SNe II, which explode first in a cold and dense medium, we allow for the 
cooling as described before. This results into an efficiency of energy tranfer of 
no more than 3% per SNII. On the other hand for type la SNe, exploding only 
after at least 30-40 Myr, we assume an efficiency of energy transfer of 100%. 

Among all the models we run, here we will present the results of the best 
ones: 
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Figure 1. The comparison between different cooling times as a func- 
tion of ISM density and metallicity. The new cooling time of Cioffi 
& Shull (1991) is shown by the solid line. The new cooling times at 
fixed Z are shown with dashed and dotted lines. The old cooling time 
indipendent of Z is indicated by the long dashed & dotted line. 



Model MG : best model of Matteucci & Gibson (1995). It is a one-zone 
model with Arimoto & Yoshii (1987, AY) IMF, cooling time and SNR evolution 
by Cox (1972) and no morphological evolution of SO galaxies allowed. We run 
it for comparison with our models with new energetic prescriptions. 

Model I : one-zone model with Salpeter (1955) IMF, Cioffi & Shull (1991) 
cooling time, SNIa without cooling and morphological evolution of SO galaxies 
allowed. 

Model II : multi-zone model (see Martinelli et al. 2000 for details) with 
the same prescriptions of model I. 

Despite of the hypothesis on type la SNR, the mean efficiency per SN 
(Ia+II) explosion is no more than ~ 20% for models I and II while it is ~ 1.7% 
for model MG. Another difference is that, thanks to the new prescriptions, the 
energy provided by SNIa makes galactic winds continuous out to the present 
time. Therefore models I and II release larger masses of elements and energy 
into the ICM than in the MG case. 



3. Time evolution of Abundances and Energy 

We consider galaxies in the range 10 9 — 2-lO 12 M0 in order to find the parameters 
linking the luminous mass Mi(z) of the galaxy and the masses of the chemical 
elements as well as the thermal energy Eth(z) ejected into the ICM. In particular, 
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we obtain relations of th type: 



E th (z) = A(z)M l m (z), 
M el (z) = B(z)M^\z), 



(3) 



(4) 



where z is the redshift and A,B , [3 and 5 are least-square fit parameters. For this 
choice of cosmological paramaters f2 m =0.3, Q\=0.7, Hq=70 km s _1 Mpc _1 , we 
integrate relations (3) and (4) over the K-band Luminosity Function in clusters 
(LF), taking into account the LF evolution with redshift (while Martinelli et 
al. 2000 did not) and the possibility of morphological evolution for SO galaxies 
into spirals for z > 0.4, for different cluster richness (n*). In order to do this 
we take the LF at z=0 in the B-band from the observational data of Sandage 
et al. (1985), then we use the transformation from B-band to K-band by Fioc 
& Rocca-Volmerange (1999) and the evolutionary corrections from Poggianti 
(1997). We consider M/L K ~ 1 at z=0 (e.g. Mobasher et al. 1999), and its 
evolution in time is calculated by means of the model of Jimenez et al. (1998). 
In order to compare model predictions with observed data we transform the 
cluster richness into temperature using the relation kT oc (n*) 8_1 . The choice 
of the K luminosity is due to the fact that it does not vary dramatically with 
galaxy evolution as it is the case for B luminosity, which is very sensitive to 
young and massive stars. 

4. Results and Conclusions 

In Figure 2 we show the iron abundance (relative to the solar value) predicted 
by model II compared to data by White (2000). An AY IMF (solid line) can 
lead to a larger amount of Fe in the ICM than the observed value. We would 
have seen the same trend if we had used model I, since one and multi-zone 
models give quite similar results (as shown also by Table 1). Our models are in 
agreement also with [a/Fe], but we predict a slighty lower O abundance. As 
we can see from Table 1, SNe can provide ~ 0.2-0.3 keV per particle. From 
the chemical point of view, we show in Figure 3 the evolution of [O/Fe] as a 
function of redshift and metallicity. At high redshift there is a fast decrease of 
O abundance due to the large production of Fe by SNIa, whereas there is no 
evidence of evolution for z < 1, in agreement with observations. 
In summary, our conclusions are: 

• In both one and multi-zone models, only those with Salpeter IMF repro- 
duce the ICM Fe abundance as well as the observed [Si/Fe] ratios. 

• Best models can provide 0.2-0.3 keV per ICM particle with a mean SN 
efficiency of ~ 20%. 

• No relevant evolution is found for abundances and energy per particle from 
z=0 out to z=l, in agreement with observational data (e.g. Matsumoto et 
al. 2000). 



• In both old and new models type la SNe play a fundamental role in pro- 
viding energy (~ 80 - 95%) and Fe (~ 45 - 80%) to the ICM. 
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Figure 2. Fe abundances in the ICM as predicted by multizone model 
compared to the observed one by White (2000), as a function of cluster 
temperature. 
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Figure 3. 



Evolution of [O/Fe] versus redshift in Coma- like cluster. 
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Table 1. Abundances in the ICM and energy per ICM particle for 
best models with evolution in time for Coma and Virgo clusters. 



Coma [O/Fe] obs. 



[Si/Fe] obs. -p- obs 



< 0.01 ±0.14 >* 



0.51± 0.60** 
< 0.14 ±0.10 >* 



< 0.23 > f 
0.33±0.05t 



[O/Fe] 



[Si/Fe] 



Fe 



Ay± gas 

10 13 M G 







E 

(keV) 



MG 
I 

II 



0.09 

-0.38 

-0.66 



0.21 

0.003 

-0.08 



-0.31 
-0.24 
-0.39 



1.27 
0.32 
0.30 



0.19 
0.13 
0.22 



Virgo [O/Fe] obs. 



[Si/Fe] obs. fs- obs. 



< 0.01 ± 0.14 >* 0.16±0.18** 
< 0.14 ±0.10 



0.40±0.02t 
0.55±0.04t 



[O/Fe] 



[Si/Fe] 



Fe 
Fe e 



M e3 
gas 

10 13 M, 



Epp 
(keV) 



MG 
I 

II 



0.08 

-0.38 

-0.66 



0.21 

0.003 

-0.08 



-0.50 
-0.35 
-0.61 



0.51 
0.13 
0.12 



0.30 
0.21 
0.34 



f (Fc/H)observcd in Coma cluster from De Grandi & Molendi (2001, Beppo-Sax);| by 
Matsumoto et al. (2000). 1(Fe/H)observed in Virgo cluster from White (2000, ASCA). 
**[Si/Fe] from Fukazawa et al. (1998). * [Si/Fe], [O/Fe] weighted mean from Ishimaru & 
Arimoto (1997) for cluster A496, A1060, A2199 and AWM7. M^ as - 1 - 20% M ICM . 
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• In spite of the large contribution from SNIa, SNe in general seem not to be 
able to provide the requested 1 keV per particle. We need perhaps other 
energy sources (e.g. quasars), although we might have underestimated the 
energy contribution from type II SNe. 

• Chandra and XMM observations on abundances and abundance gradients 
in the ICM will set stronger constraints to our models. 
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